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ABSTRACT 



Context. The structure in density and temperature of protoplanetary disks surrounding low-mass stars is not yet well known. The 
protoplanetary disks mid-planes are expected to be very cold and thus depleted in molecules in gas phase, especially CO. Recent ob- 
servations of molecules at very low apparent temperature (~ 6 K) challenge this current picture of the protoplanetary disk structures. 
Aims. We aim at constraining the physical conditions, and in particular the gas-phase CO abundance in the mid-plane of protoplane- 
tary disks. 

Methods. The light molecule HoD^ is a tracer of cold and CO-depleted environment. It is therefore a good candidate to explore the 
disks mid-planes. We performed a deep search for H^D^ in the two well-known disks surrounding TW Hya and DM Tau using the 
APEX and JCMT telescopes. The analysis of the observations are done with DISKFIT, a radiative transfer code dedicated to disks. In 
addition, we used a chemical model describing deuterium chemistry to infer the implications of our observations on the level of CO 
depletion and on the ionization rate in the disks mid-plane. 

Results. The ortho-H2D^ (li o - li.i) line at 372 GHz was not detected. Although our limit is three times better than previous obser- 
vations, comparison with the chemical modeling indicates that it remains insufficient to put valuable constraints on the CO abundance 
in the disk mid-plane. 

Conclusions. Even with ALMA, the detection of H2D+ may not be straightforward, and H^D^ may not be a sufficiently sensitive 
tracer of the protoplanetary disks mid-plane. 

Key words. Stars: circumstellar matter - planetary systems: protoplanetary disks - individual: DM Tau, TW Hya - Radio-lines: stars 
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1. Introduction 

Planetary systems are believed to be formed in the latest stages 
of protostellar evolution, within protoplanetary disks. The de- 
tails of this process are however not known. Constraining the 
physical and chemical structure in disks during their evolution 
are requirements to get more insights into this question. 

Current models predict that protoplanetary disks consist of 
three typical layers, which are, going from the surface to the 
mid-plane of the disk: 1) the outer layer, a PDR directly illumi- 
nated by the stellar UV, 2) a warm molecular zone, where the UV 
is sufficiently attenuated to allow formation of molecules, and 3) 
a cold mid-plane where temperature is l ow and most mole cules 
are expected to be frozen on dust grains (iBergin et al.ll2007l) . All 
static chemical models published so far produce similar results 



* Based on observations carried out with the Atacama Pathfinder 
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laboration between the Max-Planck-Institut fiir Radioastronomie, the 
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The JCMT is operated by the Joint Astronomy Centre on behalf of 
the Science and Technology Facilities Council of the United Kingdom, 
the Netherlands Organisation for Scientific Research, and the National 
Research Council of Canada. 
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Academia Sinica, P.O. Box 23-141, Taipei 106, Taiwan, ROC, e-mail: 
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(e.g. lAikawa & Nomural200"6tlSemenov et alJ2005h . The chem- 
istry is dominated by two main processes: molecular freeze-out 
onto grains in the cold mid-plane and photodissociation in the 
upper layers (enhanced by grain growth that allows the UV field 
to penetrate deeper in the disk). This simple picture is however 
not yet constrained by observations, and some puzzles still re- 
main. In particular, very cold C2H, CN and CO (<15 K, i.e. be- 
low t he evaporatio n temperature) were detecte d in the DM Tau 
disk (iDartois et ai]i2003: P ietu et al 20Q7; Henning et al.ll2010l: 



Chapillo n et al.l submitted) , pointing out that some process may 
lead to cold gas-phase molecules in the disk mid-plane. Radial 
and vertical mixing is usually invoked to explain partial molec- 
ular repleni shment of the cold disk m id-plane from the higher 
warm layer jSemenov et alj2006tlAikawa.2007.) . Such an expla- 
nation is not fully satisfactor y with respect to the current chem- 
ical results (see for example lHersant et al.l 12009 !). Observations 
of a chemical tracer of cold regions are thus needed. 

HaD^ is a very promising molecule to trace the disk mid- 
plane as it is likely the only remaining observable species in 
cold medium, where molecules are frozen onto grains. H2D^ is 
exclusively formed in the gas phase at low temperatures (T< 20- 
30K). Moreover, CO (and other heavy species like N2) is a 
very efficient destroyer of H3 and H2D^. High abundances of 
H2D^ can thus only build up in depleted environments. The 
ortho-H2D^ (1 1,0 - 1 1,1) line, hereafter o-H2D^, is commonly de- 
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Fig. 1. o-H2D^ spectrum in DMTau. The dashed line indicate 
the systemic velocity. The resolution is O.lVkms"' 
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Fig. 2. o-H2D^ spectrum in TW Hya. The dashed line indicate 
the systemic velocity. The resolution is 0.20 km s"'. 



tected toward s cold and dense pre-stellar cores, see for instance 
ICaselU et al. (|2008) . Obse rvation of 0-H2D+ was also proposed 
by ICeccareUi et aL (|2004|) as a way to measure the mid-plane 
ionization, since H2D^ may be the dominant ion in that region. 

Until now, no firm detection of o-H2D^ in protoplanetary 
disks has been definitively secured. A tentati ve detection on 
DMTa u obtained with the CSO was claimed bv ICeccarelli et al 
(|2004l) and later challenged by SMA observations bv loi et al 
('2008). Using a line profile appropriate for the DM Tau disk, 
[Ouilloteau et al. (2006) showed that the Ceccarelli et al. (2004) 
data indicates only a 2< t tentative detection. There are only 
upper-limits in TW Hya d Ceccarelli etani2004t iThi et aLllIOOi: 
lOi et al.ll2008h obtained with the JCMT and the SMA. 

The paper is organized as follows: the observations are pre- 
sented in Sect.|2] and upper limits on the o-H2D^ column densi- 
ties are derived in Sect.[3] The implications on the level of CO 
depletion and ionization in the disk midplane are discussed in 
Sect.|4] The feasibility of ALMA observations are discussed in 
Sect.|5]and we summarize our study in Sect.|6] 



2. Observations 

We searched for o-H2D^ in two well-studied, close or large pro- 
toplanetary disks. Their properties are summarized in Table[T] 

TW Hya is the nearest observable T-Tauri star, located only 
at 56 pc from the Sun, and surrounded by a rather small disk of 
about 200 AU of radius. It is the only protoplanetary disk where 
molecules have been studied in some systematic way because of 
its proximity (Oi et al. 2008). The disk is viewed nearly face-on. 

DM Tau is a well known TTauri st ar of 0.5 Mp surrounded by 
a large a nd massive mole c ular disk (iGuilloteau & Dutrevl 1994 
[1998; O berg et al.1 120T0I) . iDartois et all (l2003h and iPietu etM 
(1200 7) have performed high resolution multi-transition, multi- 
isotope studies of CO, and demonstrated the existence of a ver- 
tical temperature gradient in the disk as well as the presence of 
cold gas-phase CO. 

We retrieved observations of o-H2D^ toward DM Tau from 
the JCMT archive. The o-H2D+line at 372.421GHz was ob- 
served together with LI 544 in January 2007 with the multi-beam 
receiver HARP and the ACSIS backend with a resolution of 
0.03 km s '. The source was centered in the pixel at the track- 
ing center (HIO). DM Tau was observed during about 10.75 h on 
source, with a system temperature of 605 K on average leading 
to a rms of 44 mK (Tp at a resolution of 0.03 km s"' . 

Using the APEX telescope, we observed the o-H2D^ at 
372.421 GHz toward TW Hya during 3.8 h on source in July and 
September 2010 with the new FLASH 345 receiver (the dual- 
polarization receive r FLASH operates a two SB SIS mixer pro- 
vided by IRAM, see lMaier et al.ll2005h . The July and September 



observations were carried out with different backends, lead- 
ing respectively to a spectral resolution of 0.15 km s"' and 
0.06 km s"'. Data were resampled to a spectral resolution of 
0.20 km s"'. As no planet was visible, the focus was checked 
by line pointing on RAFGL 5254 and IRAS 07454-71 12, point- 
ing was done toward the same calibrators. The rms noise is 
18mK(Tp obtained after removing a constant baseline to the 
data at 0.2 km s"' of resolution. 

We did not detect any o-H2D^ emission in any of the sources. 
Assuming a Jansky-to-Kelvin (S/T^) conversion factor of 300 
for the JCMT and 41 for APE^fl and presuming that the total 
linewidths are respectively of 2 and Ikms ' for DMTau and 
TWHya, we derived 1 cr noise level on the integrated line of 
0.32 Jykms ' for DMTau and 0.33 Jykms ' for TWHya. We 
converted to flux density as th e sources are unr e solved . 

Previous observations by ICeccarelli et al.l ('2004). T hi et al.l 
( 2004) , and Oi et al. ( 20081) give a 1 cr rms of 1.2 Jy kms"> for 
TW Hya and 0.88 Jy km s~' for DM Tau (assuming a S/Tmb con- 
version factor of 55 Jy/K for the CSO telescope). So, our limits 
on the integrated intensities are respectively 3.3 and 3.6 times 
better for DM Tau and TW Hya. 

3. Results 

To derive upper limits on the o-H2D^ column density from 
the non-detections, taking into account the disk structure, we 
used DISKFIT, a radiative transfer code optimized for disks 
(Pietu et al. 2007; Pavlyuchenkov et al. 2007). This code fits a 
parametric model to the data and the best model is determined 
by a minimization process. Disks can be reasonably described 
by a parametric model where al l primary quantitie s are power 
laws functions of the radius (see iDutrev et al.ll 19941 iPietu et aH 
2007): surface density 2(r) = So(r//?o)~^, temperature T{r) = 
Toir/Royi, velocity v{r) = Voir/Ro)-" and scale height i/(r) = 
Hoir/Roy'^ that controls the vertical density structure. The tem- 
perature is a simple power law of the radius, no vertical gradient 
is taken into account. In the peculiar case of non detection, all 
parameters except the value of Eq are fixed. The surface den- 
sity of o-H2D^ at the reference radius (noted Eo) is then derived 
und er the LTE hy pothesis (valid for densities > lO^cm^^* see 
Pari se et aDl201 iL Fig. 5 a). This parametric approach limits the 
number of assumptions on the disk structure. Table|2] presents 
the resulting o-H2D^ column densities distributions for several 
models. 

For the two sources, parameters such as position, inclina- 
tion, and posit ion an gle are fixed according to previous CO ob- 
servations ([Pietu et alJ 120071: iHughes et alJ l2008l) . Two values 



' http: //docs.jach.hawaii.edu/JCMT/HET/GUIDE/het_guide 
" http: //www.apex-telescope.org/telescope/efficiency 
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Table 1. Sources coordinates (J2000) and properties 



Source 


RA 


DEC 


V,,, (kms-') 


Distance (pc) 


Radius* (AU) 


A,," (km s-') 


M. (Mo) 


Mrf,,,A-(10-^Mo)"* 


DMTau 


04:33:48.733 


18:10:09.890 


6.1 


140 


800 


2 


0.5 


3 


TWHya 


11:01:51.875 


-34:42:17.155 


2.7 


56 


200 


1 


0.5 


3 



Notes . From CO observations. FWHM derived from CO observations '***' From continuum observations dDutrev et al.lll997l ; IWilner et al.l 
I2000h assuming a dust mass absorption coeffficient k(v) = 0.1 x (v/lO'^Hz) ' cm^/g. The uncertainties are at least a factor ±3 due to the unknown 
dust characteristics. 



Table 2. o-HiD^column densities 





fixed parameters 




3(T upper limit 


Tioo 


(K) q 


Ro,„(AU) 


P 


Eioo (cm"-) 






DMTau 






15 


0.4 


720 





1.9 X 10'2 


15 


0.4 


550 


-1 


8.0 X 10" 


30 


0.6 


720 





1.1 X 10'^ 


30 


0.6 


550 


-1 


4.5 X 10" 


10 


0.0 


720 





1.3 X 10'2 






TWHya 






30 


0.5 


200 





1.4 X 10'^ 


30 


0.5 


200 


-1 


1.0 X 10'^ 


40 


0.2 


200 





1.3 X 10'2 


40 


0.2 


200 


-1 


9.0 X 10" 



Notes. Derived o-H2D^ colum n densities. Parameters similar to 
those of ICeccarelli et alj 

are adopted for the exponent of the column density distribution 
p = Q (i.e. flat distribution) and p - -\ (i.e. surface density 
linearly increa sing with radiu s) to match the results of chemical 
modeling (see IWillacvll2007l and this work). 

Given the uncertainties on the disks thermal structure, we 
derived the o-H2D^ column density under two cases for each 
source, with low (case L) and high (case H) temperature. For 
DMTau, a useful upper limit can be derived from the '^CO 
measurement with Tioo = 30K and q - 0.6. This is an upper 
limit as the '-CO emission is optically thick enough to trace 
mainly the disk uppers layers. Less abundant isotopologues, 
which sample deeper in the di sk, as well as o t her molecules 



(HCO+, CH, CN and HCN, see 


Pietu et al.l2007l;lHenning et al.1 


l201(]HChapillon et alJIsubmitted 


) give lower values (Tioo ~10 - 



15, ^ ~0 - 0.4). The law 15 X (r/100) " is adopted for the cold 
case. For TWHya, the two temperature cases considered here 
(7"ioo = 40K, q - 0.2 and Tioo - 30K, q - . 5) cor respond to 
the two temperature solutions of [ Hughes et al.l(l2008l) . 



The results are also dependent on the external radius. In the 
p - Q case, Rom is chosen in agreement with CO observations 
(i.e. 720 AU for DM Tau and 200 AU for TW Hya). When the o- 
H2D"^ column density increases with the radius {p - -1), the 
flux is dominated by the outer disk, we adopted a radius of 
550 AU for DMTau in that case. This assumption is in accor- 
dance with the "tapered edge" model (see Discussion). 

The derived 3 cr upper limits on the o-H2D^ column density 
at 100 AU is about lO'^cm^^ (ggg Table|2]). That is ~ 4 times 
lo wer than the pre v ious estimation of 4.5 + 0.9 x 10'^ cm"^ 
by ICeccarelli et al.l (l2004b on DM Tau. To make a more sensi- 
ble comparison, the valu e of Sion was also derived adopting the 
same disk parameters as ICeccarelli et al.l (2004) (Tioo = lOK, 
q=0 and p=0). In that case we find Sioo = 1.3 x 10'^ cm"^, so a 
fac tor 3 low er 

lOi et all ((2008) have derived a 3o- upper limit of 5.1 x 10'^ 
cm"^ on TW Hya. Our upper limit is a factor ~ 5 lower but be- 



cause their adopted thermal structure is not given, we cannot fit 
the same disk model to our data. 

4. Discussion 

Our deep integrations show that o-HaD^ emission, if any, is ac- 
tually much fainter than previously thought. In this section, we 
discuss the implications of our non-detections on the amount of 
CO remaining in gas-phase in the disks mid-plane. 

To derive a lower limit on the CO abundance in the disk mid- 
plane, we make use of the time-dependent chemical code dedi- 
cated to deuterium chemistry presented in Parise et al. (2011) 
and calculate the o-HaD^ abundance on a grid of densities and 
temperatures, over several hypothesis on the local CO abun- 
dance, grain growth and ionization rate. Using a realistic density 
and temperature profile for the disk, we then compute the result- 
ing column density of o-HaD^ , and compare it to the one derived 
from observations. We also derive the emission lines profiles for 
each of those disk chemical models using the radiative transfer 
code DISKFIT under LTE. 

The chemical model was first developed to desc ribe the deu- 
terium chemistry at work in cold prestellar cores (iParise et al.l 
2011). This code takes into account all different spin states of 
the different molecules and ions. In particular, ortho and para 
H2 , H3 and isotopologues are conside r ed sep arately, with up- 
to-date reaction rates from 'Hugo et al.' ('2009') and recombina- 
tion rates for all ions of Pagani et al. (2009). A simple chem- 
istry of CO and N2 is also included. The reactions are all in 
gas-phase except the formation of H2, HD and D2 on grains. 
The code is time-dependent, and the developement of deuterium 
fractionation was shown to be limited by the ortho-H2 to para- 
H2 conversion timescale (Pagani et al. 2009). In the mid-plane 
of protoplanetary disks, we expect that this conversion has al- 
ready efficiently happened, and therefore use the steady-state 
abundances reached at long times. 

We obtained the o-H2D^ abundance relative to H2 as a func- 
tion of density and temperature. We ran the code for several val- 
ues of the CO abundance (x[CO] = 0, 10"^ 10 ^). We also 
varied the effective grain size (a = 0.1, 1 and lOyum). The abun- 
dance of H2D^ should also depend on the ionization, so we var- 
ied the cosmic rays ionization rate (^ = 1 x 10"'^, 3 x 10"'^ and 
1 X 10 '^s '). The UV field is neglected, which should not be a 
problem as H2D^ will in any case be predominantly in the dense 
regions. 

To estimate the o-H2D^ column densities resulting from the 
chemical simulations we need a model of the disk H2 density 
and temperature. For the surface density we adopted two differ- 
ent models: 

- the simple model where the H2 surface density is described 
by a power law S^-(r) = S^-(r//?o) truncated at a given 
radius {Rom), hereafter Model 1, 

- and a model where the H2 surface density is an exponentially 
tapered power law as would result from viscous spreading 
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Table 3. Disk structure adopted for the chemical modeling 





H2 de 

Model 1 

E|;Mcm-2) p RaiMJ) i?„„,(AU) 


isity 

Model 2 

i:f{cm-'^) p Ro(AV) RAAV) 


Temp 
Case L 
Tm (K) q 


erature 

CaseH 
Tm (K) q 


DMTau 
TWHya 


8.4 X 10'' 0.85 45 700 
4.3 X 10-3 1.0 45 200 


9.6 X 10-' 0.45 45 180 
3.3 X 10" 0.7 45 30 


15 0.4 
30 0.5 


30 0.63 
40 0.2 



Grain 0.1 /an Grain 1 /un Grain 10 /un 

- CR 3.E-17 " 




10 1 00 1 000 

Rodius (AU) 



i 












CR 3.E- 


7 

w R 
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r 
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IT 







2 4 6 _j 8 10 

Velocity (km.s ) 



Fig. 3. Results from the chemical modeling for DM Tau, case L, 
Model 2 and ^ = 3 X 10 '^s '. Top: Predicted 0-H2D+ column 
densities as a function of radius for several values of grain size: 
a=0.1//m (black), 1 (red) and 10 (blue), and CO abundance: 
x[CO]= 10"'' (thin lines) 10"^^ (medium) and 10"^ (thick). The 
green lines correspond to the best 3cr upper limits derived from 
observations with p=0 and p=-l. Bottom: Observed spectrum 
(green) and predicted line profiles (same color and thickness 
convention). 

Table 4. Predicted 0-H2D+ line fluxes (mJykms ') for DMTau 



0-H2D+ fluxes for g = 3 x 10~"s~' 0-H2D+ fluxes for x[CO] = 10"' 



x[CO] 


0.1 


a (yum) 
1 


10 


xlO"" 


0.1 


a (jum) 
1 


10 


10"** 


14 


8 


7 


10 


3830 


932 


373 


10-5 


465 


174 


99 


3 


1730 


544 


224 


10"' 


1730 


544 


224 


10 


713 


273 


126 



of flie disk E^^r) = S^£r/^-fexp(-(r//?c)(2-p)) hereafter 
Model 2 (see also lHuehes et al.ll2008l) . 

We adopted the same temperature law as previously, without ver- 
tical gradient. Parameters are summarized in Table|3] To derive 
local densities, we further assumed hydrostatic equilibrium. 

For DMTau, the disk parameters m ainly come from the 
analysis performed by iPietu et aP (l2007h . The density models 



DM Tau 

CO abundance Cosmic ray grain size 



10' 



h ; 

□ 


; 6 

□ 




S E 

■ 

A 

- ■ \ ^ 

X 

^ , 1 


A 

■ 

A 

r ■ - 

X 

^ , ^ ^ 


m 

A 
ft 

X 

-J 1 fi-- 



W " 10 ^ 10 ' 1 _ _ lo'lO ' 1 10' 

X[CO] f (»10 s ) (ixm) 



Fig. 4. o-H2D^ predicted flux as a function of CO abundance 
(left), cosmic ray (middle) and grain size (right) for DM Tau. 
Grain size: a=0. 1 pm (black), 1 (red) and 10 (blue). Cosmic ray ioniza- 
tion rate:^ = IxlO"" (triangle), 3x10"" (square) and lxlO""'s"' (oc- 
tagon). CO abundance: x[CO]= 10"' (open marker) 10"^ (filled marker) 
and lO"'' (cross). The green line corresponds to our 3cr limit derived 
from observations. The dashed lines illustrate the dependencies. 

are der ived from the continuum observations of lGuilloteau et al.] 
Model 1 is arbitrarily extended up to an external radius 
of 700 AU (in agreement with CO emission). 

Fo r TW Hya, the parameters are derived from iHughes et al.l 
(l2008h . Note that their "cold" case applies to the density Model 
1 and the "hot" case to the Model 2. 

The disks are expected to have vertical temperature gra- 
dients which we ignore in this simple modeling. As a conse- 
quence, the abundance of o- H?D^ above the mid-plane i s under- 
estimated (see for instance lAsensio Ramos et al.ll2007l Fig. 1). 
Nevertheless, as the H2 density above the mid-plane is signifi- 
cantly lower, the contribution of this additional layer to the total 
0-H2D+ column density is negligible (see Fig. 6 and 7 of lWillacvl 
2007). For DMTau the two adopted temperature profiles repre- 
sent two extreme cases, one only considering the low mid-plane 
temperature, the other using the higher values sampled by '^CO 
at 1-2 scale heights. 

The top panel of Fig. [3] presents the distribution of the o- 
H2D^ column density for DM Tau, Model 2 case L and the stan- 
dard value of the cosmic ray ionization ^ = 3 x 10"'^ s"'. 
All models with CO depletion show a region of high column 
density at R = 40-50 AU (low temperature case) or R = 150- 
200 AU (high temperature case). In each case, this corresponds 
to a kinetic temperature around 20 K. These models display 
also a region of low o-H2D^ column density for T~15K. This 
is caused by the rapid evolution of the ortho/para ratio of H2D^ 
with temperature in this range (the ratio o/p decreases by a 
factor of 3 between 20K and 15K for a density of 10'' cm"^, 
and x(CO)=10 "''). This eff'ect has a lready been sho wn in previ- 
ous w orks, see iFlower et al.l (|2004 their Fig. 6) and lSipila et al.l 
(l2010i their Fig. 4). Apart from that feature (which is much 
more pronounced for x[CO]=0), the o-H2D^ column density in- 
creases ■mtlrthe radius in Model 1 (in qualitative agreement with 
IWillacvlbOOTh . In Model 2 this increase is tapered off by the 
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exponential fall of the H2 surface density. For DM Tau, the o- 
H2D^column density peaks around R - 550 AU where it reaches 
id'^cm -2. The disk of TWHya being small (R„„, ~ 200 AU), 
the temperature is high in the whole disk, so the o-H2D^ column 
density is low whatever the CO abundance (column density be- 
low lO'^cm"^ if we exclude the high ionization case) and the 
models are not constraining. The predicted integrated spectra for 
each of the models are computed using DISKFIT. The bottom 
panel of Fig. [3] presents the predicted spectra together with ob- 
servational data for DM Tau, Model 2 case L and the standard 
value of the cosmic ray ionization. The predicted o-H2D^ line 
fluxes for DM Tau are presented in Table |4] and Fig|4] The in- 
tegrated line flux seems to vary roughly like and a""-^. The 
dependancy over the CO abundance is more complicated. 

From Fig.|4]one can see that unless the ionization rate is very 
high, the current upper limits are not very stringent on the disk 
parameters. At best only the small grains, CO depleted by a fac- 
tor 100 and ^ > 3 x 10 '^s ' case are excluded in the DM Tau 
case. With a CO depletion of order 10 or so, any grain size would 
fit the current data. The chemical models are not constraining for 
TWHya. 

Co mpared to the earlier study of ICeccarelli & Domini^ 
(l200l . which did not consider the ortho and para H2D^ sep- 
arately, our chemical model including all spin-dependent re- 
action rates predicts column densities on average 5-10 times 
weaker. Accordingly, as the o-H2D^ line emission is mostly op- 
tically thin, except for the highest column densities (above a few 
10^^ cm"^ ), w e predict flux densities 3 to 5 times lower than 
lAsensio Ramos et al.. C2007.) . 



5. Implications for ALMA 

In this section, we discuss the feasibility of observational im- 
provements compared to the present work, using the ALMA in- 
terferometer. Figure 3 demonstrates that an increase in flux sen- 
sitivity by at the very least a factor of 10 (note that the emission is 
optically thin) compared to the present observations is required 
to secure a convincing detection of H2D^ (in the framework of 
the chemical models of Sect 4). An even higher sensitivity would 
be required to disentangle the different models. The current up- 
per limits, obtained with single-dishes of 12 - 15 m diameter and 
integration times 5-10 hours, show that only larger instruments 
can expect detecting H2D^. ALMA is the most promising: in 
single-dish mode, the 50 12-m antennas will bring a gain of 7 in 
sensitivity compared to APEX (assuming similar receiver per- 
formances). However, in interferometric mode, even in its com- 
pact configuration, the angular resolution of the main array will 
be around 1", so that the sources will be highly resolved out. In 
such cases, the brightness sensitivity becomes the relevant pa- 
rameter In the beam of the current instruments, the observed 
discs are diluted by about a factor 3 to 4, as the APEX beam 
size is 14" compared to an outer disc diameter around 7 - 8". 
However, the aperture filling factor of ALMA is only about 40 
%, compared to filled aperture for APEX. Hence, in a similar 
integration time, and assuming similar noise performance, the 
ALMA brightness sensitivity should enable to reach a detection 
level of column densities only 2 times smaller that the current 
level in each ALMA synthesized beam. Further spatial smooth- 
ing would improve the sensitivity to extended structure, but only 
moderately, as this implies degrading the weight of the longest 
baselines. Using the ALMA compact array (ACA) may be more 
appropriate: the collecting area of the 12 7-m antennas is 4 times 
that of APEX, so accounting for an effective filling factor around 



30-40 % also leads to similar brightness sensitivities as for the 
12-m antennas compact configuration. 

While detailed predictions will require proper knowledge of 
the ALMA performances at the difficult frequency of 372.4 GHz, 
these simple considerations show that o-H2D^ is unlikely to be- 
come a major "workhorse" to study e.g. the kinematics of disks. 

6. Conclusion 

- We have improved the sensitivity to o-H2D^ by a factor 3 
over previous observations. No o-H2D^ is detected, in con- 
tradiction with previous tentative detection. 

- With our chemical model, TW Hya is in any case too warm 
or not dense enough to provide detectable o-H2D^ at the cur- 
rent sensitivity. 

- The limit on DM Tau is incompatible with very high CO 
depletion (x[CO]< 10"^) and small grains (0.1 mum). 
However, it is not discriminant with more reasonable values 
of the grain sizes and/or CO abundances. 

- Much more sensitive observations will be required to obtain 
significant constraints. Even with the powerful ALMA inter- 
ferometer, it will be extremely difficult to reach a sensitiv- 
ity that will allow to constrain the ionization level and the 
amount of CO in the disk mid-planes by means of the obser- 
vation of o-H2D^. 
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